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LONGSHORE SAND TRANSPORT STUDY AT CHANNEL 
ISLANDS HARBOR, CALIFORNIA 


by 
R.O. Bruno, R.G. Dean, C.G. Gable, and T.L. Walton, Jr. 


I. INTRODUCTION 


1. Purpose. 


The relationship between longshore sand transport in the surf zone and wave 
energy is of vital interest to coastal engineers concerned with design and main- 
tenance of navigation and beach erosion control projects. As discussed in the 
Shore Protection Manual (SPM) (U.S. Army, Corps of Engineers, Coastal Engineering 
Research Center, 1977), past field and laboratory studies have produced a widely 
used empirical relationship for sand transport in a shore-parallel direction; 
however, these studies were conducted in areas where total transport may not 
have been measured. In this study, Channel Islands Harbor, California (Fig. 1), 
was selected for measuring longshore sand transport because an offshore break- 
water and twin jetties at the site form a unique sand trap. This site is 
considered to be nearly a total littoral barrier to longshore sand transport 
moving in a southerly direction (Herron and Harris, 1966). A further advantage 
to the Channel Islands site is its exposure to high wave energy climate with a 
dominant wave direction out of the west and northwest and consequently, high 
sand transport rates in a predominantly southerly direction along the beaches. 
Dredging records from the harbor show annual transport in excess of 1 million 
cubic meters (Herron and Harris, 1966). The objective of this study was to 
obtain additional data on longshore sand transport in order to reevaluate the 
present semiempirical relationship for longshore sand transport given in Chapter 
4 of the SPM. 


2. Location and History. 


Channel Islands Harbor is located about 110 kilometers northwest of Los 
Angeles and 1.6 kilometers upcoast (northwest) of Port Hueneme (Fig. 1). The 
harbor consists of entrance jetties and an offshore breakwater which were 
constructed to provide a solution to a downcoast beach erosion problem and to 
provide a small-craft harbor. Before the harbor construction, a serious 
erosion problem existed downcoast (southeast) of Port Hueneme which has been 
attributed to diversion of littoral sands into the Hueneme Canyon by the Port 
Hueneme north jetty. Presently, sand trapped at Channel Islands Harbor is 
bypassed to the south of Port Hueneme by periodic dredging (U.S. Army Engineer 
District, Los Angeles, 1948; Herron and Harris, 1966). 


The offshore breakwater is 700 meters long, is located 600 meters offshore 
in a water depth of 9 meters, and trends roughly parallel to the shore. The 
design of the sand trap was developed empirically by considering the configu- 
ration of the Santa Monica, California, breakwater fillet and by developing 
diffraction patterns for generally prevailing waves (U.S. Army Engineer 
District, Los Angeles, 1948; Herron and Harris, 1966). 


36 Physiography and Littoral Processes. 


The shoreline at the site forms the coastal edge of the Oxnard Plain, an 
abandoned flood plain of the Santa Clara River (Fig. 1). The low, flat flood 
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Figure 1. Channel Islands Harbor study area. 


plain extends about 21 kilometers along the shoreline and 13 kilometers inland, 
and is bounded on the north by the Sulphur Mountains and on the south by the 
Santa Monica Mountains. These mountains terminate at the sea in hard, wave- 
resistant formations, forming the south bank of the Ventura River and Point 
Mugu, respectively. The principal drainage features are the Ventura and Santa 
Clara Rivers. Offshore slopes in this area are gentle except where the steep- 
walled Hueneme and Mugu submarine canyons cut the Continental Shelf to within 
0.4 kilometer of the shore. 


Figure 2 is a view of the study area in relation to offshore islands and 
shows that much of the wave action is intercepted or modified by Point 
Conception or the offshore islands. The principal avenues of wave approach 
are from the same direction, hence, sea and swell arrive predominantly from 
the northwest and west. Local winter storms of short duration and a limited 
amount of summer swell, originating from the south Pacific Ocean, reach the 
Hueneme area from the southwest and create short periods of northward littoral 
drift. However, wave studies and long observation of the shoreline processes 
show conclusively that there is a great preponderance of southward littoral 
drift (Herron and Harris, 1966). 


_ There is one source area of littoral material, in addition to updrift 
beaches, for the Channel Islands study site (Fig. 1)--the Santa Clara River 
which discharges at the upper end of the Oxnard Plain. At irregular intervals 
of 10 to 30 years, tremendous floodflows occur that form a large delta at the 
mouth of the river extending as much as 0.8 kilometer seaward of the normal 
alinement of the shore (Herron and Harris, 1966). However, there is no accurate 
measure of the rate at which littoral material is supplied by the river. 


At Port Hueneme, the Hueneme submarine canyon extends to within about 300 
meters of the shore and is the one major sink for littoral materials in the 
Channel Islands area. The profile from mean lower low water (MLLW) to the 
-18-meter depth steepens from a normal of 1 on 100 to 1 on 4 into the canyon. 
The steep slopes of this canyon continue to depths as great as 1.5 kilometers. 


Past estimates of littoral sand transport have shown a net southward sand 
transport on the order of 920,000 cubic meters per year for the beach area 
between the Hueneme and Mugu submarine canyons, based on shoreline erosion 
surveys before the Channel Islands bypassing operation (Herron and Harris, 1966). 


The entrance to Channel Islands Harbor was dredged in 1960-61, and since 
that time the material impounded in the sediment trap has averaged approximately 
1 million cubic meters per year; the capacity of the trap is approximately 2 
million cubic meters. The material is usually pumped (on a biennial basis) a 
distance of about 2 kilometers to the downdrift side of the entrance to Port 
Hueneme, although material is occasionally placed in the "pocket beach" formed 
by the southeastern jetty to Channel Islands Harbor and the northwestern jetty 
to the Port Hueneme entrance approximately 1.6 kilometers to the southeast 
(Herron and Harris, 1966). 


The net sediment transport is from northwest to southeast, but periods of 
reversal do occur. The average sediment size at the site is approximately 
0.20 millimeter, the textural characteristics of the material deposited in the 
trap were reported in Bruno, Watts, and Gable (1977). 
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One feature of an offshore breakwater is that it tends to "overtrap" 
material, especially when waves approach in a direction nearly perpendicular 
to the shoreline. This effect can be demonstrated by noting that even for waves 
propagating directly toward shore, the effects of diffraction will cause depo- 
sition in the form of a salient or tombolo behind the breakwater. If this 
deposited material is removed such that an equilibrium topography does not 
exists, more material will be transported behind the breakwater, etc. Unfor- 
tunately, the degree of overtrapping for different wave characteristics is 
currently unknown. A related question of whether the impounded volumes should 
be correlated with the net or downdrift component of longshore energy flux is 
discussed in Section II. 


For more detailed analyses of shoreline changes and historical bypassing 
operations in this area, see Beach Erosion Board (1953), Savage (1957), and 
Herron and Harris (1966). 


4. Past Longshore Sand Transport Relationships and Studies. 


The relationship between immersed weight transport rate of sediment, I, 
and longshore energy flux factor at breaking, Pgs, is generally presented in 
terms of a dimensionless constant, K (e.g., Vitale, 1980). 


I = KPog (1) 


where for a single wave train of breaking height, H,, and breaking direc- 
tion, Cpe relative to a normal to the beach, Pp. can be expressed as 


yf 
Pes = —g © sin ap cos ap (2) 
in which y is specific weight of water, and C the wave group speed at 
breaking. A more useful form of equation (1) is obtained by noting that the 
immersed weight transport rate, I, and bulk longshore sediment transport, 
Q, are related by 


I 


eee ee (3) 


in which p and pg are the mass densities of water and sediment, respec- 
tively, g the gravitational constant, and P the inplace porosity of the 
sediment. Combining equations (1) and (3) 


gc eyenie dant! te a1 
(== - 1)ga - p) 


The above equations for longshore sediment transport do not account for 
sediment size, beach slope, porosity, etc.; therefore, it is very unlikely 
that K is a true constant. Tests conducted to date have not provided the 
ranges of independent variables (e.g., sediment size) or accuracy necessary 
to distinguish the effects of these variables, although models which account 
for the effects of these variables have been postulated (Dean, 1973; Walton, 
1979). 


Bs = te 
Q = Pe Pe (4) 


The right side of equation (4) is the manner in which the sand transport 
relationship is presented in the SPM, where K* fis assumed constant and has 
dimensional units. 


Available field data indicate that none of the approximately 56 data points 
are based on both measured wave direction and volumetric accumulations in a 
total trap. Some are based on visually observed wave directions at a location 
quite distant from the measured wave heights. In addition, much of the sedi- 
ment data are based on tracer studies with the attendant uncertainties in 
estimating depth of effective motion. 


The following brief review of field data on sand transport relationships, 
used in the SPM and other studies (see Galvin and Schweppe, 1980), is helpful 
in assessing the currently used relationship. Laboratory data have been 
omitted as they do not represent a properly scaled version of the prototype, 
since the scaled sediment size is significantly larger than is usually present 
with sand-sized particles. The review is presented in chronological order of 
the various investigations which yielded the data. Other reviews include Das 
(1971) Greer and Madsen (1978), and Bruno, Dean, and Gable (1980). 


a. South Lake Worth Inlet, Florida (Watts, 1953). Longshore transport 
rates were based on quantities of sand transferred by a permanent sand by- 
passing plant on the north jetty of South Lake Worth Inlet. The pressure drop 
of the bypassing pump was correlated with sand discharge by a series of pumping 
events into a diked disposal area on the downdrift (south) side of the inlet 
and subsequent surveys of the associated volumes. Thereafter, transport rates 
were inferred from the (calibrated) pressure drop of the pump. Wave charac- 
teristics were based on-(1) wave height measurements from a staff gage mounted 
on the South Lake Worth pier about 16 kilometers north of the inlet, and (2) 
visual observations of wave direction at the surfline as obtained from a van- 
tage point about 6 kilometers north of the inlet. This study yielded four 
data points with a sediment diameter of approximately 0.4 millimeter. The 
average K value is 0.90 with a standard deviation of 0.11. 


b. Anaheim Bay, California (Caldwell, 1956). Dredged material from the 


entrance to Anaheim Bay was placed on the downdrift (southeast) shore, and 
repeated surveys of this area were conducted as the material was transported 

in a southerly direction. Changes in volume were interpreted as longshore 
transport rates and estimates of longshore component of wave energy flux were 
based on (a) wave staff measurements from the Huntington Beach pier about 9 
kilometers to the south, and (b) wave directions based on hindcasts and recog- 
nition of the sheltering by the offshore islands for waves originating from 
certain directions. This study provided five data points with a sediment diam- 
eter of approximately 0.40 millimeter. The average K value is 0.76 with a 
standard deviation 0.38. 


c. Cape Thompson, Alaska (Moore and Cole, 1960). The growth of a spit 


and associated waves were observed over a 3-hour period. Spit volumes were 
measured by plane table survey, and wave characteristics were based on visual 
estimates. Only one data point was obtained with a K value of 0.25 for a 
sediment diameter of 1.00 millimeter. 


d. Silver Strand, California, and El Moreno, Baja California (Komar and 
iE 


nman, 1970). These data represent transport over fairly short-time intervals 
as determined from sand tracer measurements and wave energy fluxes derived 


from an array of wave sensors. Sand transport volumes were usually determined 
over a fraction of a tidal cycle as the product of the width of the surf zone, 
the longshore displacement of the center of gravity of the tracer, and the 
thickness of tracer movement. The latter quantity was based on observations 
of the depth to which a cylindrical "plug" of tracer had eroded over the obser- 
vational period and depth of tracer in cores; depth values ranged between 2 
and 10.5 centimeters. These measurements yielded a total of 14 data points, 
and the average and standard deviation values of K for Silver Strand (4 data 
points) are 0.77 and 0.18, respectively; the average sediment size was 0.18 
millimeter. These values for El Moreno (10 data points) are 0.82 and 0.27, 
respectively; the average sediment size was 0.60 millimeter. 


e. Channel Islands Harbor, California (Bruno and Gable, 1976). Data in 
this report were based on the same general field program that is the subject 
of the present report. Sediment transport rates were inferred from volumetric 
accumulations behind the Channel Islands Harbor offshore breakwater, and the 
longshore energy flux values were based on Littoral Environmental Observations 
(LEO). A total of 13 data points resulted with average K and standard 
deviation values of 1.61 and 1.19, respectively. Sediment size was approxi- 
mately 0.2 millimeter. These data were not used in the SPM sand transport 


relationship. 
5. Discussion. 


The results presented above comprise a total of 37 data points which, 
excluding Bruno and Gable's (1976) data, are summarized in Table 1 and 
presented in Figure 3. Figure 3 represents the SPM sand transport rela- 
tionship in the dimensionless form of the equation. Figure 4 attempts to 
discern any effect of sediment size D on the quantity kK. There is a 
fairly reasonable relationship of increasing K with decreasing sediment 
size, although there is only one data point for a sediment size exceeding 
0.6 millimeter. 

In general, it appears that none of the data sets provide the confidence 
that should be associated with both the wave and sediment volume data. In those 
studies with well-established volumetric data, the wave data usually included 
one or more visually estimated wave parameters (height and direction). In the 
reported tracer studies, it is believed that the inferred sediment transport 
rates overestimate the actual transport. This expected bias is probably due 
to the estimate of the shore-parallel displacement of the center of gravity of 
tracer displacement being obtained from the upper 5 centimeters of the sediment 
column, since shear transport must exist within the sand bed with the upper 
layers moving most rapidly. Additionally, the use of a single value of trans- 
port thickness based on the erosion depth of a tracer plug should increase bias 
as the maximum erosion depth is expected to increase with time. Finally, the 
longshore sediment transport in the surf zone is expected to be both spatially 
and temporally variable although these scales are currently unknown. Since the 
tracer studies represent longshore sediment transport over only a part of a 
tidal cycle, it is surprising that the studies do not exhibit greater scatter, 
i.e., standard deviations of 23 and 33 percent for the Silver Strand and El 
Moreno K values, respectively. One advantage of a relatively long-term, com- 
plete sediment trap is that it integrates the temporal and spatial variability, 
thereby providing a good basis for investigating the mean structure of the trans- 
port phenomenon before attempting to measure and understand the fine structure. 
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Figure 3. Comparison of field data with present 
SPM relationship. 


© Band G, N=!3 


\ 
\ 
\ 
\ 
\ 
\ 
\ 
\ 
\ 
\ 
N 
N 
\ 
W,N=4e. 
N ® Kand!, N=I0 
@ e \ ? 
Kand!|,N=4 \ 


Figure 4. 


0.5 1.0 
t, D (mm) 
Plot of K versus diameter, D. Initials and numbers, N, 
refer to investigators and number of data points comprising kK. 


I 


The complexity of longshore sediment transport, the uncertainties and dif- 
ferences exhibited in available data, and the economic and functional impact 
of a valid quantitative predictor of transport certainly justify substantial 


future field and laboratory investigations. 
II. DATA COLLECTION AND ANALYSIS 


The data collection program consisted of periodic bathymetric and topo-— 
graphic surveys, routine wave data from which longshore transport and the wave 
energy flux could be calculated, and sediment samples taken during the study i 
period to obtain quantitative information on the sand size at the site. Much 
of the information used here was reported in Bruno and Gable (1976). 


The wave data included both LEO (Berg, 1968; Bruno and Hiipakka, 1972) and 
gage data from two wave pressure sensors installed at the site. A chart de- 
scribing the various types of data collected and the survey dates from 1974 to 


1977 is presented as Figure 5. 
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Figure 5. Longshore sand transport study, Channel Islands, 1976-77. 


1. Bathymetric Survey Data Collection and Analysis. 


U.S. Army Engineer District, Los Angeles assisted in the surveys. A base 
line was established approximately parallel to the shore, and profile lines 
normal to the base line were spaced at about 30.5 meters for a distance of about 
823 meters (Fig. 6). Elevations at each station were measured from the base 
line out to the detached breakwater. Surveys were scheduled at 4- to 6-week 
intervals, but the scheduling was frequently modified due to survey crew avail- 
ability, equipment failures, and unfavorable surf conditions. 


Table 2 provides a summary of survey data, showing survey dates, ranges 
surveyed, type of fathometer calibration, survey method, and qualitative esti- 
mate of overall data quality, Under survey method, "standard" indicates an 
analog fathometer record was made, and the Los Angeles District's standard sur- 
vey techniques were used to position the survey vessel. A man onshore directed 
the vessel operator along the profile line. At 10-second intervals, the analog 
recording was marked and the vessel position was recorded by plane table and 
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Figure 6. Longshore transport study area. 


Table 2. Summary of survey data. 
Ranges surveyed Fathometer Data 

calibration quality 

101+00 to 122+00 Leadline Standard 
101+00 to 127+00 Leadline Standard 
101+00 to 127+00 Leadline Standard 
101+00 to 127+90 Leadline Standard 
101+00 to 127+00 Leadline Standard 
101+00 to 127+00 Leadline Standard 
101+00 to 127+00 Leadline Standard 
101+00 to 127+00 Leadline Standard 
101+00 to 122+00 Leadline Standard 
101+00 to 127+00 Leadline Standard 
101+00 to 127+00 Bar check Hybrid 

101+00 to 127+00 Bar check Hybrid 

101+00 to 127+00 Bar check Hybrid 

101+00 to 127+00 check Standard 
101+00 to 127+00 check Hybrid 

101+00 to 127+00 check Hybrid 

101+00 to 134+00 check | Automated 
101+00 to 134+00 check | Automated 
101+00 to 134+00 check Automated 
101+00 to 134+00 check Automated 
101+00 to 134+00 check Automated 
101+00 to 134+00 check Automated 
101+00 to 134+00 check Automated 
101+00 to 134+00 check Automated 
101400 to 134+00 check Automated 
101+00 to 134+00 ae check Automated 
101+00 to 134+00 check Automated 
101+00 to 134+00 check Automated 
101+00 to 134+00 check Automated 
101+00 to 134+00 check _ Automated 
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The "hybrid" method indicates an analog fathometer record was made and 
the survey vessel positioned by use of electronic ranging equipment. With this 
method, a man onshore directed the vessel operator along profile lines, but the 
vessel position was monitored by telemetering data to a field office where a 
real time plot of position was produced by a minicomputer. Monitoring the sur- 
vey in progress determined and eliminated both errors in positioning. An elec- 
tronic timer at 2-second intervals simultaneously marked the analog record and 
recorded the vessel's position on magnetic tape. 


The "automated" method indicates that the soundings from the fathometer 
were digitized electronically and transmitted with position information to the 
field office minicomputer. All survey data soundings as well as position were 
monitored in real time. 


A LARC 5 amphibious vehicle was used in all three methods for the bathy- 
metric part of the survey. The LARC allows the continued measurement of the 
survey line through the surf zone, using rod and level methods when the 
vehicle's wheels contact the bottom. Profile lines over the dry beach areas 
were measured by standard level and rod transects. 


Late in phase I of the study, it was found the standard method was not 
providing the reliable, accurate data desired. As shown in the plots in Figure 
7, the data are on a line sheltered by the breakwater and at a point beyond 
normal sand deposition. Divers measured the bathymetry using underwater ref- 
erence stakes to verify no deposition at this point. These plots indicate that 
unacceptable errors were introduced on several surveys. These errors were a 
result of poor positioning data, poor fathometer calibrations, and errors in 
data reduction. 


The standard surveying method assumes that the LARC is on line for all fixes 
and that the vessel's position and the fathometer record have been marked both 
precisely and simultaneously. However, this is not always true and errors that 
occur can never be recovered. Even with no errors produced in field collection, 
reduction of working scale plane-table data sheets was limited to an accuracy of 
about 3 meters which was seldom attained. 


In the hybrid method, an electronic timer simultaneously marked the fathom- 
eter at 2-second intervals and controlled the recording of the LARC's position, 
eliminating human errors in determining position. Position accuracy, which was 
limited by the electronic ranging equipment, was about 2 meters and consist- 
ently attained. 


During the surveys, errors were also produced in fathometer calibrations 
and data reduction procedures. Under the standard method, leadline soundings 
were used for the first year (April 1974 to March 1975) for fathometer calibra- 
tion. In taking leadline soundings, accuracy is determined by numerous factors, 
including skill and care of the operator, motion of the vessel, and roughness 
of the water's surface. Later surveys used a bar lowered at 5-foot intervals 
under the fathometer for calibration. This type of calibration (bar check) 
eliminates human factors and is complete in that it shows any occurrence of 
nonlinearity in the fathometer. 


Uncertainty in interpreting the fathometer records is another source of 
error which was discovered. Figure 8 shows fathometer records of two profiles, 
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both with similar features which appear to be waves; however, on one record 
these features indicate the irregular bottom left after dredging. For this 
reason, all data reduction under the hybrid method was performed by the field 
investigators who observed conditions at the time of data collection. Also, 
this possible misinterpretation was the reason no smoothing was done under the 
automated method. 


The data quality rating in Table 2 is a subjective rating, from 0 (poor) 
to 10 (good), determined by the field investigators after examining the data on 
hand, field notes, and all profile plots. Using this approach, it was deter- 
mined that data collected 18 June 1974 and 27 March 1975 (see Fig. 7) were 
unreliable (Table 2). 


In the survey analysis, deposition from upcoast (northerly) transport 
during times of reversals in wave direction and possible deposition due to 
reversals in the area of the trap seaward and upcoast of the jetties were 
examined. Figure 9 is a plot of the first profile upcoast of the jetties. 
Survey dates of 17 April 1974 and 16 September 1975 are plotted, showing the 
total change at this profile over that period. The plot shows little deposi- 
tion at the end of the jetty. All profile data were examined in detail to 
define the area of deposition attributable to influx of material from north 
(upcoast) to the trap. No significant deposition was measured at distances 


more than 485 meters from the base line. Using this distance (485 meters) from — 


the base line as an outward boundary to calculate volumes, errors due to influx 
of material from reversals were minimal. 


Volume changes were calculated, using the area between the base line and 
485 meters offshore and the north jetty and 670 meters upcoast. Over this area 
a volume accuracy of within +5,000 cubic meters was expected when survey errors 
had been minimized. 
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Figure 9. Range 101+00 at beginning and end of study. 
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Figure 10 shows bottom changes at station 109+00 (midprofile of the trap 
area) between April 1974 (immediately after dredging in phase I of the study) 
and July 1974 (4 months later) surveys. Predredging conditions in September 
1973 are also shown. 
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Figure 10. Profile changes at range 109+00. 


Figure 11 shows bathymetric changes within the sand-trap area from 17 April 
to 30 July 1974; only major changes in bathymetry are reflected due to con- 
touring at l-meter intervals. The map indicates a zone of deposition from sta- 
tions 102+00 and 122+00 along the shore and 30 to 300 meters seaward of the 
working base line. The figure is representative of the type of filling the 
sand trap experienced during the entire project study. Most of the sand-trap 
deposition occurred well within 300 meters of the shoreline. 
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Figure 11. Elevation changes between 17 April and 30 July 1974 surveys. 
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2. LEO Data Collection and Analysis. 


Figure 6 shows three sites where twice daily surf data were collected using 
procedures developed under the CERC LEO program (Bruno and Hiipakka, 1973; Berg, 
1968). These data include observations of surf conditions, local winds, and 
littoral currents. The primary concern of this study was the estimates of 
breaker height, breaker direction, and longshore current. To aid in estimating 
the breaker direction, LEO observers were provided with a protractor on the 
data form. Longshore currents were measured using small packets of dye which 
disperse upon injection in the surf zone. The observers measure the distance 
the dye travels parallel to the shoreline. Current speed is estimated from the 
movement of the dye patch centroid over a 1-minute period, and current direc- 
tion is noted. Surf zone current velocities are not uniform; therefore, the 
width of the surf zone is estimated as well as the distance from the shoreline 
to the point of dye injection. 


LEO data were used to compute the longshore energy flux factor, using two 
methods of computation, in order to relate the wave energy to the longshore 
sediment transport. The first method involved the typical approach to computing 
the longshore energy flux, using wave heights and wave angles at breaking. This 
method utilizes the SPM equation (4-29): 


pglp® : 
Pima SS 16. © sin 2ap (5) 
where 

Pos = longshore energy flux factor 
f°) = fluid density 
Hy = breaking wave height 
Op = breaking wave angle 
C = group wave velocity at breaking. 


Using linear wave theory, shallow-water assumption and a linear relationship 
between breaking wave heights and depth equal to 


Pecan) (6) 


dy 


where dj} = breaking wave depth, equation (5) becomes 


5/2 WD 
Ps = ob : (2) sin 2a (7) 
The longshore energy flux computed from equation (7) is dependent on only two 
surf zone parameters--breaking wave height and breaking wave angle. The 
assumed value of « for computation, k = 0.78, is reasonable for various mono- 
chromatic wave theories (e.g., see Longuet-Higgins, 1970). The breaking wave 
height measured by LEO observers is assumed to be the significant wave height 
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(average of the highest one-third of the waves). In this study it was desired 
to compute the "true'’ wave energy flux; therefore, wave height was modified to 
reflect the root-mean-square (rms) wave height in accordance with the statis-— 
tical wave train theory of Longuet-Higgins (1967) assuming a narrow-banded wave 
energy spectrum. For this particular situation, which results in a Rayleigh 
probability distribution for wave height, the relationship between significant 
wave height and rms wave height is as follows: 


Hepa) (2) auriene (8) 


where H, is the significant wave height, and H;m, the root-mean-square wave 
height. 


The second method computes longshore energy flux, using breaker wave height 
and longshore current as determined from LEO measurements, along with estimated 
values of surf zone width, W, and distance from shoreline to injection point 
of dye, X. The equation for determining longshore energy flux using this 
second method is 


pg Hy WV Ce 


Pos = 5t K vi (9) 
2 alle 


where V is the longshore velocity as measured in LEO program, Cr¢ the friction 
coefficient dependent on bottom water particle excursion bottom roughness, and 

(V/Vo)LH the theoretical dimensionless longshore velocity for Longuet-Higgins' 

assumed mixing parameter, P = 0.4; (V/V,);y is dependent on the parameter X/W 

(see Longuet-Higgins, 1970). 


The derivation of equation (9) is given in Appendix A. The rms wave height 
is used in equation (9) to calculate the longshore energy flux. Longshore 
velocity is typically measured over a timespan on the order of the "modulated" 
wave train period. 


In this part of the study, additional emphasis was placed on determining a 
reasonable bottom friction coefficient, C¢, for use in equation (9). The fric- 
tion coefficient was determined by minimizing the variance between computed 
longshore currents (from LEO observation) and theoretical longshore currents 
(from Longuet-Higgins' (1970) formula with an assumed mixing coefficient 
P = 0.4). LEO data from the two northernmost LEO stations (out of the zone of 
breakwater influence) were chosen for the calculation. More than 2 years of 
twice daily wave observations (4,464 observations) from the two stations was 
used to compute the friction coefficient, Cs = 0.0056, which was used in equa- 
tion (9) to calculate the longshore energy flux. 


In both methods, weighted averages of the longshore energy flux at breaking 
are calculated which correspond to periods between volumetric surveys. Values 
of Pgs are given for the three stations in Appendix B, Tables B-l and B-2. In 
most survey periods, the value of Po, computed from current data was higher 
than the value of Pgg calculated using wave angle information. In addition 
to longshore energy flux, correlation coefficients R are also calculated for 
the correspondence between calculated longshore currents (from Longuet-Higgins' 
(1970) theoretical approach) and measured longshore currents. Correlation 


coefficients ranged from -0.18 to 0.76 and are independent of assumed friction 
factor. 
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In assessing the data for use in predicting sand transport, it is pertinent 
to assess the quality of the LEO data from the standpoint of consistency, both 
spatially and temporally. Values of Pp, computed from LEO stations 5714 and 
5715 (Fig. 6, two northernmost stations) show surprisingly good agreement when 
comparing similarly computed values within the same time frame (either current 
or wave angle). This correspondence also holds for station 5713 which is only 
90 meters north of the north end of the Channel Islands offshore breakwater. 
It is believed that station 5713 observations were strongly affected by dif- 
fraction effects caused by the breakwater. Station 5713 data were not used in 
any of the correlations between the longshore energy flux and sand transport 
discussed later in this report. 


To indicate temporal consistency of the data source, average Py. values 
(for stations 5714 and 5715) were plotted versus time (in months) for the two 
phases of the study (Fig. 12). Comparing the two phases in the figure shows 
that the low wave energy periods occur consistently during the summer and early 
fall months. 
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Figure 12. Net longshore energy flux versus time for LEO 
stations 5714 and 5715. 


Correlation between the weighted average longshore energy flux and the long- 
shore sand transport is discussed in a later section of this report. Problems 
arose when reversals in wave direction (and consequently the longshore energy 
flux) occurred, due to the limited extent of the survey area in which accre- 
tion was always dominant. As a result, there were periods in which both the 
weighted average longshore energy flux and the volume change were positive 
(northward-directed energy flux and accretion). This apparent discrepancy was 
due to the limited extent of the survey area north of the offshore breakwater 
which was still trapping sand even though northward of the survey area, the 
sand transport was in a northerly direction. This overtrapping effect is due 
to wave diffraction patterns caused by the breakwater driving sand into the 
sheltered area behind the breakwater. In these particular cases, no correlation 
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could be made. Thus, the effects of a strong reversal in wave climate must 
also be considered in the other data since an evenly balanced wave climate with 
a preponderance of southward-directed energy flux (which would give a low nega- 
tive "net" P,. value) might be correlated with an excessively large volumetric 
accretion. To assess the possibility of this phenomena, the ratio of the cumu- 
lative total northward energy flux to the cumulative total southward energy 
flux was calculated for the average Py, value of stations 5714 and 5715. This 
ratio is presented in Appendix B, Table B-3, and is discussed later in this 


report. 


At the end of phase I of the study, LEO wave observations were compared to 
the wave gage significant wave heights to determine any inconsistencies in the 
LEO recorded waves. Figure 13 provides the comparison, for a 6-month period 
during phase I of the study, which shows that the observer consistently re- 
corded lower wave heights than the gage recorded in the case of the higher 
waves. Although the plotted gage wave heights have not been modified to re- 
flect a breaking wave height, an analysis of the transformation to breaking 
conditions in a few selected cases did not alter the results significantly. 
This observer bias in underestimating the higher waves was found to result from 
the observer practice of recording the last breaking wave closest to the shore- 
line in the case of multiple-breaking waves. Thus, when the higher waves broke 
offshore, re-formed, and then broke again, the observer was only recording the 
last and closest (to the observer) breaking wave height. 
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Figure 13. Comparison of observed wave heights (LEO) to recorded 
wave gage heights in phase I of study. 
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In phase II of the study the observer had been requested in estimating wave 
heights to record the outer breaker. As a result, a comparison of LEO wave 
height measurements and wave gage measurements in phase II gave a much better 
correlation (see Fig. 14). 


3.0 
LEGEND 


LEO within | hr of Gage Measurement (82 Obsns.) 
LEO within |-2 hr of Gage Measurement (106 Obsns.) 
LEO within 2-3 hr of Gage Measurement ( 95 Obsns.) 


2.5 


283 Observations taken during Dec. 76 - Aug.77 
aot Channel Islands Harbor, California 


2.0 
E 
= 
isp 
3 
iB a3) 
> 
tS) 
= 
no} Att 
@ 
fe 
2 4 + 
Co) I} (0) yee Bieta 
® ® 6 nX.;) 6A + 
A eg Sakis 3 a 
c og sarah, +h eh, SECO 
eo + @ ce} 
+ @ + A @A vA 
@ A+ one @deded + +0 
+@ ® 
7 
0.5 7 8 


7 +h4 A ® ® 


4 9@ 


0 0.5 1.0 1.5 2.0 2.9 
Significant Gage Height (m) 


Figure 14. Comparison of observed wave heights (LEO) to recorded 
wave gage heights in phase II of study. 


The strong dependence of the longshore energy flux on wave height 
Geile = H5/2 | eq. 7; Pyg ~ H, eq. 9) necessitated disregarding the LEO data in 
phase I and using only the wave gage directional spectra for sand transport 
correlations. 


In comparing the calculated values of longshore energy flux from the two 
methods, the wave angle method (eq. 7), and the longshore current method (eq. 
9), it was found that the longshore current method led to higher values of 
longshore energy flux during most data periods. The reason for this is unknown. 
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3. Wave Gage Data Collection and Analysis. 


During the period 30 July 1974 to 6 May 1975, two total pressure sensors 
were located near the bottom in an approximate 6-meter water depth, separated 
by approximately 23 meters; the geometric characteristics of the gage instal- 
lations are presented in Figure 15. Pressure records, which were recorded in 
digital form approximately every 2 hours, contained 4,096 data points at 0.25 
seconds each, resulting in a record length of 17.07 minutes. In order to cal- 
culate Pe, at breaking, the following steps were carried out: (a) Calculation 
of the frequency-by-frequency wave direction and energy at the location of the 
wave gages; (b) transformation of the wave spectrum to the breaker line, in- 
cluding shoaling and refraction effects; and (c) computation of Py, at the 
surfline. Each of these is described below. 


y 


x (normal to contours) 
Gage 1 
Figure 15. Characteristics of a two-gage array at Channel 
Islands Harbor, California. 


a. Calculation of Wave Characteristics at Wave Gages. The pressure time 
series were analyzed using a standard fast Fourier transform (FFT) program to 
determine the coefficients. For example, the time series from gage 1 is rep- 
resented by 


ee i2tnj 
‘ [2,@ - ib, (mle % (10) 


P,(t;) = S; 
n= 


in which i = v-l1, and N is the total number of data points (N = 4,096) for 
each gage. 


The FFT coefficients are defined in terms of the pressure time series as 

i2tnj 

ays , 

aj (n) - ib, (n) = N 2 Py(t;)e (11) 
ik 


and the a,(0) term represents the mean of the record. In calculating the FFT 
coefficients, there are a large number of "tapers" (or windows) that may be 
used to reduce the adverse effects of spectral leakage which arises due to 
representing an aperiodic time series by a periodic series (Harris, 1976). 
These tapers all have the form 


(Zs) 


(Gea) = w(j)p(t,) (12) 


in which p(t;) is the measured pressure, and w(j) is a weighting factor. A 
characteristic of these tapers is that they are unity at the midpoint of the 
time series and decrease to a lesser value near the two ends. In the present 
analysis, it was determined through a number of test calculations that although 
the effect of tapering at an individual frequency could be substantial, the 
overall effect on Pgg was quite small, generally less than 4 percent. Thus, no 
taper was used in the analysis. 


The depth of water, Ad, overlying the pressure sensors is obtained from 
a, (0) and ay(0) as 


ad = 9 [ay (0) + a9(0)] (13) 


in which y is the specific weight of seawater. The total water depth, d, 
is the sum of Ad and the distance, Sc, of the pressure sensors above the 
bottom (0.3 meter). 


Each FFT pressure coefficient is transformed to a water surface displace- 
ment coefficient by the following linear wave theory relationship: 


[an(n), bn(n) } = TO) [ap (n), bp(n)] (14) 


in which kK, @) is 


cosh k SQ) 


cosh kd (5) 
n 


K, (n) = 


In equation (15), ky, is the wave number associated with the angular frequency, 
On, as obtained from the linear wave theory dispersion relationship 


0, = gk, tanh k,d (16) 


One disadvantage of measuring waves with near-bottom pressure sensors is 
evident by examining equations (14) and (15). For higher frequencies (shorter 
wave periods) Kp is very small, resulting in very small pressure fluctuations 
near the sea floor for the higher frequency waves. Thus, to avoid contaminating 
the calculated water surface displacement, n, it is usually necessary to apply 
a high-frequency cutoff above which the pressure contributions are discarded. 
The proper selection of this high-frequency cutoff depends on the signal-to- 
noise characteristics of the pressure sensor and signal conditioning system. 

For the present analysis, the high-frequency cutoff was established at a wave 
period of 3.1 seconds. For a nominal 6-meter water depth and a 0.3-meter height 
of the pressure sensor above the bottom, the pressure signal is attenuated to 
approximately 16 percent of its surface value. In the survey area, where rea- 
sonably long Pacific swells occur, neglecting wave energy for periods shorter 
than approximately 3 seconds appears justified. ; 
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Denoting hereafter the FFT coefficients for the water surface displacement 
as a(n) and b(n), it is noted that the coefficients have the following prop- 


erties: 
Sau SNS 
iS 
n= 
a 
Bos ie 
b 
N 
5 +k 
and thus 
heh N/2 
qe a 2 oF 
n=1 


Therefore, the total (kinetic 
a particular frequency component, 


E(n) = 


[a2(n) + b?(n)] Gly) 
ay 
- $-% 
(18) 
-by 
= aia 
[a2(n) + b?(n)] (19) 


plus potential) energy, E(n), associated with 
ny els 


2y[a2(n) = b*(n)] (20) 


For wave direction, consider the definition sketch in Figure 16 and the 
following representation for n(x,y,t3): 


N SE ee cos a(n) cos a(n) + ky (n) sin a(n) - e(n)] 
n(x,y,t;) = » Fie 
n=1 
. a 
SR ([eG) © alyGwy ile — & (21) 
n=1 
y t 
Gage 2 / 
J 4 
/ 
Sy 
/ 
4 
Mi Alternote Wove Crest 


Wove Crest 


wave pirection 


a (shoreword ) 


Figure 16. The two-gage array notation and directional ambiguity. 
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Denoting the water surface displacements at the two gage locations as 
n(x,,y,,t) and n(xo,y2,t) and calculating the cross spectrum, Si2(n), 


S}2(n) = [a,(n) - iby(n)] [a2@) - ib2(n)] (22) 


where the overbar denotes the complex conjugate, 


Syo(n) = [az (n)az(n) + bj (m)b2(n)] - ilag(m)by (nm) - a; (n)b2()] 


= c19(n) - iqi2(n) (23) 


and cy2(n) and q,2(n) are the "cospectrum" and "quadspectrum," respectively of 
ny and nz. S}9(n) can also be expressed as 


-i[k, (xo = x}) "eos a(n) 4 kj, Go = yq)) sin c@))] (24) 
Sy)o(n) = F*(n)e 


and noting from Figure 16 that the separation distance and orientation of a 
line joining the two gages are denoted by & and 8, respectively, S,7(n) 
can be written as 


=ik fzcos [a(n)a=78) 


Sradn) = B>Gne 
(25) 
= F2(n) cos {knt cos [a(n) - 81] - iF2(n) sin{kn2 cos [a(n) - e1| 
Comparing equations (23) and (25), it is apparent that the wave direction, 
a(n), relative to the x-axis can be expressed as 
& Sah Gym GY) 
a(n) = B + cos ! ean 1) Receee (26) 
ky,e C15 (n) 


where the + is the result of a directional ambiguity associated with a two- 
gage array (Fig. 16) and is avoided by choosing the minus sign which selects 

the wave arriving from the seaward half-plane adjacent to a line joining the 

two gages. 


There are two conditions for which it was not possible to calculate the 
wave direction, a(n): (1) poorly conditioned wave data, presumably due to 
spectral leakage, and (2) spatial aliasing due to the fairly large separation 
distance (~23 meters) between the two gages. If the first condition exists, 
the absolute value of the quantity within the brackets {} in equation (26) may 
exceed unity, clearly a physically impossible condition since the extreme values 
of the cosine function are +1. This tended to occur for very long waves for 
which both the energy and the value of ka were small. The percentage of 
energy within an individual record when this condition occurred was relatively 
small, averaging 2 to 3 percent with a maximum of approximately 10 percent. 
The second condition requires that the wavelength be equal to or greater than 
twice the projection of the wave gage separation distance in the direction of 
wave propagation. Referring to Figure 16 
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L > 22[cos(a - B) de as (27) 


which indicates that for the least adverse effects of spatial aliasing, the 
gages should be alined parallel to the dominant orientation of the wave crests. 
In the present case 6 = 86.8° and, in general, the waves approached within 
approximately 40° of the shoreline. In the analysis, a variable aliasing fre- 
quency limit depending on the wave directions for the lower frequencies was 
developed for each record. Three frequency limits were considered, corre- 
sponding to la - B| values of 45°, 60°, and 75°, with the 75° value associated 
with the highest frequency limit. Selection of the frequency limit required 
that 90 percent of the wave energy in a frequency range O,y/j to oq lay with- 
in the associated wave direction range, where o, is the frequency corre- ' 
sponding to the particular aliasing wave direction (eq. 27). For the water 
depth at the site, the three directional limits of 45°, 60°, and 75° are associ- 
ated with approximate wave periods of 5.3, 4.2, and 3.1 seconds, respectively. 
The percentage of energy associated with periods shorter than 5.3 seconds was 
substantial in the data, amounting to 25 to 30 percent. In later calculations 
of Pp,, an attempt was made to account for the energy above the aliasing fre- 
quency by augmenting the calculated values, using 


ETOT 


Vgc SE Bieuers (28) 
cm (© 


in which the subscripts c and cm indicate calculated and calculated modified, 
respectively. Epo; and E represent the total wave energy values and the 
energy below the spatial aliasing frequency. This modification is equivalent 
to associating the effective direction as determined from the frequencies not 


affected by the aliasing consideration or poorly conditioned data to all the 
wave energy. 


b. Determination of Breaking Depth. At this stage, the wave energy and 


wave direction in the vicinity of the gages are determined. These values are 
then transformed to the breaker line, accounting for shoaling and refraction. 


To determine a breaking depth, a method was used in which the total onshore 
flux of wave energy was equated at the gages and at the breaker line with the 
requirement that the rms breaking wave height, Hrms, be related to the water 
depth by 


Hrms S Kdy (29) 


in which k was taken as the usual spilling breaker value of 0.78. The total 
onshore energy flux, Fp, at the reference (gage) location is 


in = 2p 2 fae@) & De @y Coin) 6°88 %n (30) 


(n) 


in which the subscript R denotes the reference location. Considering shallow- 
water conditions for breaking Cg (n) = C= ved. and cos a, ~ 1.0, the onshore 
energy flux at breaking is 


33 


2 0 
Ved, = Ye Jed (31) 


Equating Fp and Fp) 


16 5 [a2(m) + b?(n) IC y(n) cos Op 
BA Ws ih Innis haha et A NL SS (32) 
ae 


.Note that neglecting wave refraction in the breaking depth determination does 
not result in very large errors; e.g., for a wave direction at the gages of 
30°, the error in breaking depth would be less than 6 percent. For a more 
realistic overall wave direction of 20°, the associated error in breaking depth 
is less than 3 percent. 


c. Transformation of Wave Components to Shore. With the breaking depth 


known, each wave component is transformed to shore, accounting for both wave 
refraction and shoaling as based on linear wave theory. Wave refraction was 
computed in accordance with Snell's law and the assumption that straight and 
parallel contours existed between the gage and breaking locations, 


Cy (n) 


ap (n) = re 


Cat) sin ia (33) 


Shoaling was based on linear theory, resulting in the value of the sums of 
the squared FFT coefficients at the breaker line of 


cos a, (n) [C (n)] 


2 2 Soe Sea CLEN NG 
lary) = De) Is See %(n) [Cy ()], 


[a*(n) + b*(n)]p (34) 


in which the first and second ratios on the right side of the equation repre- 
sent the effects of refraction and shoaling, respectively. 


d. Computation of Pp. at the Surfline, with the wave energy and direction 
known at the breaker line, the value of the wave energy flux factor, Py,, is 
readily determined. 


N/2 
Piet = Given [a2(n) + b*(n) Jy Cop (n) [cos a(n) sin a(n) ]} (35) 
n=1 


where the factor G is given by the ratio 


@ 6 —SS= (36) 
as defined in and discussed in relation to equation (28). 
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e. Orientation of Wave Gages. In calculating directional wave properties 
relative to the shoreline, it is extremely important to establish the orienta- 
tion of the wave gage array and a representative shoreline orientation. During 
the original study, the locations of the wave gages had been established by 
measuring angles with transits from two known shore stations on the project base 
line. These angles were then translated into gage locations through standard 
geometric procedures. Two sets of such measurements were taken on 24 and 26 
June 1974. The resulting gage characteristics are presented in Table 3. 


Table 3. Geometric characteristics of the two-gage array. 
Gage separation distance | Gage orientation with 
(m) respect to base line 


24 June 1974 23).55, 102° 
26 June 1974 22.71 14.44° 
7 June 1980 

Teneul, I 23.8 15 (Cul 

ietail D DBD LAC TES 

wei 3 22.6 1236s 

Values adopted 


in this study 


=e 


IS}, 2B? | 


Although the gage separation distances obtained from the two surveys dif- 
fered somewhat more than desired, the gage orientations that differed by 3.6° 
were of particular concern. Thus, on 7 June 1980, these variables were estab- 
lished using range-range microwave equipment. Measurements were obtained from 
three combinations of shore-based locations, which allowed a better determina-— 
tion of gage orientation; the results are also presented in Table 3. The average 
of all values were adopted for the present study. 


f£. Effective Shoreline Orientation. The determination of effective shore- 
line orientation was somewhat more difficult at this site than originally antic- 
ipated. One difficulty was that dredging and subsequent filling by longshore 
sediment transport in the impoundment area caused local anomalies in the bottom 
contours, such that they were not quite straight and parallel. In particular, 
following a dredging activity, the nearshore contours moved landward more than 
the seaward contours. This is believed to be an expression of the concentration 
of the longshore sediment transport in the nearshore zone. 


The shoreline orientation was established by plotting and overlaying pro- 
files located in the vicinity of the wave gages. The relative landward-seaward 
locations of the profiles were adjusted in the overlays to obtain best agree- 
ment by eye and the associated orientation determined. A range of orientations 
were obtained by matching the nearshore and offshore parts of the profiles. 
These ranges were qualitatively consistent with those determined from aerial 
photography, with the larger angles associated with the more seaward parts of 
the profiles (see Table 4). 
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Table 4. Shoreline orientation as determined from 
beach and offshore profiles. 
Date of survey Range of effective shoreline 
orientation relative to base line 


8 May 1974 9.95° to 11.33° 
18 June 1974 8.83° to 10.78° 


Value adopted 
in this study 


4. Results of Longshore Sand Transport-Longshore Energy Flux Factor Correlation. 


In the correlation of Pp, with I, only the values directed toward the trap 
were used. The rationale is that during all periods when waves were directed 
toward the trap, there was sediment available for transport; therefore, this 
variable provides the best basis for correlation. 


The volumetric accumulations were interpreted as I values by equation (3), 
using p = 0.35, and by dividing the time (in seconds) between surveys. 


In phase I of the study (April 1974 to September 1975), the available values 
of longshore energy flux factor were computed using wave gage data. Longshore 
energy flux factors were also computed using LEO data during the same period 
but, as discussed previously, for higher waves the LEO wave heights, in com- 
parison with wave gage data, were typically underestimated so the LEO data were 
not used in phase I. The wave gage data were available nominally every 2 hours. 
The value of Py, for a survey period was obtained by summing all Pg, values 
directed toward the trap and dividing by the total number of wave recordings. 


In phase II (January 1976 to August 1977), only data from one wave gage 
were available for a limited time; therefore,.it was not possible to analyze 
directional wave spectra to compute the longshore energy flux factor. As dis- 
cussed earlier, the LEO wave height observations during phase II did correlate 
reasonably well with wave gage data. As a result, the longshore energy flux 
factor was computed using LEO data and the LEO current method presented in 
equation (9). 


Table 5 provides a summary of the analyzed data in phases I and II. Figure 
17 presents the same data along with the SPM volumetric sand transport relation- 
ship converted to the more common relationship of equation (1) between the 
longshore energy flux factor and the immersed weight sand transport rate (in 
metric units). In this relationship, the constant K is dimensionless and 
ranges from 0.42 to 1.51 in the phase I data and from 0.19 to 4.18 in the phase 
II data. The average value for the two sets of data is K = 0.98 with a break- 
down for the two sets as follows: 


Phase No. of data K °K °K/K x 100 

iin ated ty bye axadg ead y43 aa biDEty) we 
I 7 oxe7"29"0.33 38 
II 11 LaOS) 1) sO) 103 
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Immersed Weight Transport Rate, I (N/s ) 


Figure 17. 


Table 5. Summary of analyzed data.!>? 


Survey period I Pes Pog = Wave data 
(N/s) | (N/s) | (N/s) SOUFcE 


Phase I 
30 July to 20 Aug. 1974 
20 Aug. to 24 Sept. 1974 
24 Sept. to 6 Nov. 1974 
6 Nov. 1974 to 7 Jan. 1975 
7 Jan. to 11 Feb. 1975 
11 Feb. to 4 Mar. 1975 
4 Mar. to 14 Apr. 1975 
Phase II 
20 Apr. to 20 May 1976 
20 May to 28 June 1976 
28 June to 12 Aug. 1976 
12 Aug. to 8 Sept. 1976 
8 Sept. to 6 Oct. 1976 
6 Oct. to 1 Dec. 1976 
1 Dec. 1976 to 3 Feb. 1977 
3 Feb. to 22 Mar. 1977 
22 Mar. to 6 June 1977 
6 June to 21 July 1977 


gages(2) 
Wave gages(2) 


Wave gages(2) 


Wave gages(2) 


Wave gages(2) 


Wave gages(2) 


Wave gages(2) 


LEO 


INegative values of longshore energy flux factor, P,,, repre- 
sents flux toward the sand trap. 


Pos values using wave data were calculated in accordance with 
equation (35); Py, values using LEO observations were calculated 
using equation (9). 


(SPM, 1977) 


LAT 


LEGEND 


@ Phase L Pes Values Computed Vio Wove 
“Gage Directional Spectra 


© Phose IX Py, Values Computed Vio LEO 
Current Method 


\! i 
' 2 3 4567810 20 3040 60 80100 200 400 600 
Longshore Energy Flux Factor at Breaking, Pos (N/s ) 


ship for Channel Islands study. 
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1,000 


Sediment transport energy flux relation- 


Figure 18 presents the present data set along with the data set used in the 
SPM relationship. 


1,000 
ae ee eS aa 
[cesar | ress cee] [| jes 
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mE eel en 

200 


LEGEND 


Watts, D:0.4mm 

Komar and Inman, D=: 0.18 mm 
Komar and Inman, 0:0.6 mm 
Caldwell, D- 0.4 mm 


Immersed Weight Transport Rate, 1 (N/s ) 


This Study, Pgg- Via LEO 
Current, D=0.2mm 

This Study, Pg, Via Wave Gage 
Directional Spectra, D-0.2 mm 


1 2 3 4567810 20 30 40 60 80100 200 400 600 1,000 
Longshore Energy Flux Factor at Breaking, Ps (N/s ) 


Figure 18. Summary of field data I versus Pp,, including 
results of this study. 


Using the data in the SPM sand transport relationship along with the results 
of the present study, an average K can be found, based on 42 field data 
points, to be K = 0.87. 


This value of K is only 11.5 percent higher than the existing value of 
K in the SPM. In view of the large standard deviations encountered in the 
data sets, it is suggested that the existing SPM relationship remain unchanged 
until additional quality data are analyzed in future studies. 


5. Conclusions. 


Measurements of immersed sediment transport rate, I, as inferred from 
sediment accumulation in a near-total trap were correlated with longshore energy 
flux factor at breaking, P 5? as calculated in phase I from two pressure sensors 
mounted near the bottom an in phase II from LEO wave height and longshore cur- 
rent observations. Based on the analysis of 7 correlation periods between sand 
accumulation and wave gage data (phase I) 11 correlation periods between sand 
accumulation and LEO wave observations (phase II), it is concluded that: 
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(a) On an overall basis, the results of the present study agree 
reasonably well with some of the previous studies. The dimensionless 
constant K in the relationship I = KPy, was found to range from 0.19 
to 4.18 with an average value K = 0.98 for the Channel Islands data. 
In view of the large scatter within the data set, it is recommended 
that no change be made in the SPM (U.S. Army, Corps of Engineers, 
Coastal Engineering Research Center, 1977) relationship for sand 
transport. 


(b) For the phase I survey periods, the K values obtained here 
are substantially lower than those based on visual wave estimates 
obtained by Bruno and Gable (1976). Their average K value is 1.40 
as compared to 0.98 obtained here for the Channel Islands data. 


(c) A pair of wave gages appears to perform reasonably well to 
determine overall directional characteristics such as Ppg,; however, 
the quality of the detailed directional results could be improved 
substantially by a greater number of gages. 


(d) In calculating wave direction from an array of gages, it is 
extremely important to accurately determine the relative locations 
of the gages and the shoreline orientation. For the present data, 
an angular misalinement of +0.5° can cause an associated change in 
Pes, ranging from 8 to 18 percent. 


(e) LEO wave height observations appear to be subjective with 
strong observer bias. As a result, LEO measurements should be cor- 
related with wave gage or other instrumentation measurements before 
their acceptance. 


(f) There are variations in the K values that cannot be explained 
through any evident mechanisms or processes other than poor data 
quality in specific data sets. 


(g) The environmental and economic significance of an improved 
quantitative predictor for longshore sediment transport justifies 
substantial future field programs to discern the effects of individ- 
ual variables (e.g., grain size, beach slope, etc.) and in general, 
to develop a more reliable predictor. 
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APPENDIX A 
DERIVATION OF EQUATION (9) FOR LONGSHORE ENERGY FLUX 


a. From Longuet-Higgins (1970) 


Vy = te (Ge) C@Ap)?/? sin ap cos ap) (A-1) 
where 
B = a mixing parameter (dimensionless) 
dp = breaking depth 
m = beach slope 
a = breaking wave angle 


b. Using relationship 


Cs SS (A-2) 


equation (A-1)} becomes 
_ 50 /«8\ _1(/gHy m sin 2ab) 
) v« (A-3) 


c. Longshore velocity at any point within the surf zone can be defined as 


V = Vf a Nh (A-4) 
» (NE; 

where V is the longshore current within surf zone and Vo, the theoretical 

longshore velocity at breaking with no mixing. 


d. From equation (58) in Longuet-Higgins (1970) 


Vo eh (A-5) 
Wey B 
e. Using equations (A-3), (A-4), and (A-5), longshore velocity is 
Vv \/5t\/K 1) — 
y aEale (Cr EEG, SED Achy) (A-6) 
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f. Using equation (4-40) in the SPM 


2 
og H i 
Pha = eo C sin 2ap (A-7) 
where C = group celerity = Yed. according to linear shallow-water wave theory; 
therefore, 
Hy) fide \u 2 
pg “b b 3 
Bushs siep ys (z) veH, sin 20, Ce 


g. Using equations (A-2), (A-6), and (A-8), and assuming m = W/dp, the 
following equation is found: 


ce (A-9) 


h. The value of (V/Vg) can be assumed equal to that given by Longuet-Higgins 
(1970) 


he aaa 


i. Equation (A-9) now becomes 


Hy, WVC 
b OSreabe uaket: (A-11) 
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APPENDIX B 


CALCULATIONS OF LONGSHORE ENERGY FLUX FACTORS 
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Table B-3. Ratio of cumulative longshore energy fluxes using 


LEO data. 
Comparative Dates Obsns. | Ratio=ZPp ,+/ZPp .- 
survey 


Phase I 


17 Apr. to 7 May 1974 

7 May to 30 July 1974 

30 July to 20 Aug. 1974 

20 Aug. to 24 Sept. 1974 
24 Sept. to 6 Nov. 1974 

6 Nov. 1974 to 7 Jan. 1975 
7 Jan. to 11 Feb. 1975 
11 Feb. 1975 
4 Mar. to 14 Apr. 1975 
14 Apr. to 6 May 1975 
6 May to 1 Aug. 1975 

1 Aug. to 13 Aug. 1975 
13 Aug. to 16 Sept. 1975 


to 4 Mar. 


oO ON DN F&F WW DY 


ay 
oO 


Poe 
NR 


13 


Phase II 
14 


15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 


28 Jan. to 20 Apr. 1976 

20 Apr. to 20 May 1976 

20 May to 28 June 1976 

28 June to 12 Aug. 1976 

12 Aug. to 8 Sept. 1976 

8 Sept. to 6 Oct. 1976 

6 Oct. to 1 Dec. 1976 

1 Dec. 1976 to 3 Feb. 1977 
3 Feb. 1977 
22 Mar. to 27 Apr. 1977 
27 Apr. to 6 June 1977 

6 June to 21 July 1977 

21 July to 30 Aug. 1977 


to 22 Mar. 


to the north (right when looking offshore); negative longshore energy 
flux factor signifies energy flux to the south (left when looking 
offshore). 


48 


Le9 BAG} Pe darecn: €0COL 
é Coe! 
sou $ ((°S9%n) Jequa) yoreessy SutTiseuTSuq Te3seoj) aeded [TeoTuyse] 
rSeTIIg "II *O°y founag *T *ASi1aue aAeM *y *saAem Aq Jiodsue13 
pues *¢ *sassao0id Te1073T] °Z °(°FTTeED) SpueTS] Touuey) °*] 
*qaodsueiq Juoawtpes 
aIOYSSUOT pue SOTISTIIIIEISYD BAeM UdEMJeq SUOT}JETAIIOD sUTWASIep OF 
suoTjeAtesqo TenstaA Aq pue siadnpsuerj sansseid pajunow woj}0qG-1e8U 
om} Aq apeu atom Squomainseay *eTUIOJTTe) S10qieYyY spueTS] TeuueYD 
ye Apnjs aeak-z e ut pouteqqgo (JJTIP [Te1073TT) seyoeeq BuoTe j10d 
-suei] pues jo uotqotpead 10j3 poyjem pejepdn ue saptaoad j1ioday 
wT86T TFady,, 
*aTIT} AaAg 
(Z-1g Ou § AaqUeD YoIeesey BuTIseuTsuq 
Teaseoy Amay °s*g / aoded Teotuysey)---wd /Z : *TIT : *d [gy] 
*1861 “SIIN Wozjz 
aTqetreae : °ea ‘SpteTy3utads { azaque) yoieessy ZutirseuTsuq Te IseoD 
Awdy *S*n : °e&a SATOATOEG JA04--*[°Te 3a] °° OunIrg *O°y Aq / eTU 
-1OJITH) ‘10qieY spueTS] Touuey) ae Apnys jiodsueiqz pues aioussuoyT 


109 CoC nOU darecn: €02OL 
Cols 
sou § ((*S9°f) Jequey yoreessay BuTAseuT3uUq TeISeOD) aeded TeoTUYyde] 
:SeTzseg "II *o°y Sounag *] *ASi9ua aaemM *H *saaem Aq Jaodsuelq 
pues °¢ ‘*sassad0id [er07IT] °*Z *(°FTTED) SpueTs] Teuueyp *] 
*jiodsueij Juewtpes 
aIOYSBUOT pue SOTISTAAIIeTeYD BAeM UdEMJAq SUOTIETEIIOD sUTWIZIepP OF 
suoTjeAtesqo Tensta Aq pue siaonpsuei} sansseid pejunow woz}O0qG-1e98U 
om} fq opel e10oM SjJUaWeINSeaW *eTUIOFTTe) *10qiey spueTS] TouueYD 
ye Apnjs izeak-z © ut poeutej3qo (3JTAP Te10IITT) Sayoeeq BuoTe aod 
-suei} pues jo uotTjotpeiad aoz poyjzew pajepdn ue saptaoad jiodey 
wT861 THady,, 
*aTITI APAOCD 
(Z-1g *ou § JajqUeD YoTeasey BuTIseUuTsUq 
qeyseoy Awszy seq / aaded Teotuysey)--°wo 77 : *TTt : *d [gh] 
“T1861 “STIN wo1z 
atqeTqeae : °ey ‘ppetz8utads { raque9 yoivesey BuTrseuTsuq TeqseOD 
Away *S°n : °eA SatoATeg 3aog--*[*Te 32a] °°* ounag *O°y Aq / eTU 
-1OJITe) ‘10qiey spueTSs] TeuueyD 3e Apnqs jaodsue1q pues s10yssu0cyT 


£29 7-18 “ou datecne €07OL 
2o= 18 
sou £ ((°S°N) JequeD YyoIeessy BuTAseuTSUq TeqIseo)) deded TeoTUYydIe] 
:SeTJ9S "II °O°yY founag *] *ABiaua oAeM *h *sSaAem Aq Jiodsueij 
pues *¢ ‘*sessad0id Te1072T] °Z “*(°3TTeD) SpueTSs] TouueyD °T 
*qaodsueij Juawtpas 
@IOYSZUOT pue SOTISTIeqVIeIeYD aAeM UdEMJeq SUOTJeETSeIIOD sUTWIaJap 07 
suoT]BAIaSqO TensTA Aq pue siadoNpsueaj osansseid pajunow woz.0qG-1e9EU 
oma Aq apew a1eM SqUoMeINsSea “eeTUIOJTTe) f10qIeY spueTS] TeuueYD 
qe Apnqs iaeak-zZ e ut pautej3qo (JJTAIP Te103ITT) seyoeaq BuoTe j1aod 
-sueij pues Jo uot OTpead 10j¥ poyjew pajepdn ue saptaoiad jioday 
w°T86T THady,, 
*3TITI eA 
(Z-1g ‘ou § i1equep yoIeesey BuTrseuTsuy 
qTeaseo) Away °s°n / ateded Teotuysey)--°wd /Z : *TTE : sd [gy] 
“1861 ‘SILIN wo2jz 
eTqeTTeae : ea ‘ppTetysutads { aequa yoreaessy BuTAseuTsuyq Te seo) 
Awzy *S°n : °ea S1TOATegG J1Oq¥--*[°Te 3a] °°° ounag *Q°y Aq / eTU 
-1IOJITS) ‘10qiey spueTs] Touuey) ye Apnys jiodsue1ijz pues ai0ys3u0y 


£29 7-18 *ou dargsn° €079L 
sea18) 
sou § ((*S°M) 1eqUaD YoIeessy BuTiseuTsuq Te}seo)) aeded TeoTuYyoey 
:S®TIeS “II ‘O°y Sounzg °] °AB19ue aaeM “yh *saaem Aq JIOdsue1y 
pues *¢ ‘*sasseoo0id Te103IT] *Z °(°FTTeD) SpueTs] Teuuey) °] 
*qaodsueij JueutTpes 
aAIOYSSUOT pue SOTISTAGVIOeIeYD BAeM UdEMJEq SUOTIeTIIIOD suTWAIaIep OF 
suoT}eArasqo Tensta Aq pue siaonpsuei} eansseid pajunow wozj,0qG-1e8U 
om] Aq apem atom Squolloinsea *eTUIOJTTeQ S10qiey spueTS] TeuueYy) 
qe Apnjs 3zeak-z e ut peuteqqo (3JTIP T[e1033TT) seyoereq B3uoTe qa1od 
-suei1} pues Jo uotjotpead 103 poyjzem pejepdn ue saptaoad jaoday 
.T86T THady,, 
*aTITI sc 
(Z-1g cou { Jaqueg yoieasay BuTiseuT3uq 
Teaseoy Away °s°en / aeded Teoptuysey)---wo /Z : “TIT: *d [84] 
°1861 ‘SIIN wo1z 
atqetrTeae : ce, ‘prtetsy8utads $ aaque) yoreassy BuTasouT3ug Te seo) 
Away °S*n : cea S1TOATOg 3IO¥--*[°Te Je] °°* ounag *Q°y Aq / eTU 
-1OJTTe) ‘10q1ey spueTS] Toeuuey) ae Apnjs yaodsue1q pues s10yssu0y 


ae 


ts 


